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We present a physical scenario in which both Fermi arcs and two-dimensional gapless Dirac 
states coexist as boundary modes at the same two-dimensional surface. This situation is realized in 
topological insulator-Weyl semimetal interfaces in spite of explicit time reversal symmetry breaking. 
Based on a heuristic topological index, we predict that the coexistence is allowed when (i) the 
corresponding states of the Weyl semimetal and topological insulator occur at disconnected parts of 
the Brillouin zone separated by the Weyl nodes and (ii) the time-reversal breaking vector defining the 
Weyl semimetal has no projection parallel to the domain wall. This is corroborated by a numerical 
simulation of a tight binding model. We further calculate the optical conductivity of the coexisting 
interface states, which can be used to identify them through interference experiments. 


Introduction —Protected surface states are the key 
characteristic of topological phases of matter. Time re¬ 
versal invariant topological insulators (TI) host at their 
surface two-dimensional (2D) massless Dirac (masiparti- 
cles protected by time reversal symmetry (T)P^. Weyl 
semimetals (WSM) are 3D gapless materials described 
at low energy by Weyl fermions. They host topologically 
robust surface states referred to as Fermi arcs since they 
form an open Fermi surfac^^HH^ Their emergence can be 
understood in terms of charge conservation (gauge invari¬ 
ance) of the effective field theory describing the WSM’s 
response to external electromagnetic fields, in analogy 
with the quantum Hall effecfP^HI^. Such an effective re¬ 
sponse predicts a number of striki ng physical properties, 
such as a finite Hall conductivit}®^^^^^^^^, a current 
par allel to an external magnetic field (chiral magnetic ef- 
fect)P^^^^I^^II^, and a finite angular momentum induced 
by a thermal gradient (axial magnetic effect 

The band structure of a WSM is characterized by a 
linear dispersion around a set of non-degenerate band 
touchings points called Weyl nodes. Their existence re¬ 
quires the breakjM of either time reversal T or inver¬ 
sion X symmetry^. Each Weyl node is chiral and has 
an associated momentum space Berry flux that gives it 
the character of a Berry flux monopole. Since the to¬ 
tal flux in momentum space is required to be zero by 
gauge invarianc^^, the Weyl n odes m ust appear in pairs 
with opposite monopole chargJ^^^. They are therefore 
topologically stable as they can only be annihilated by 
bringing together a pair with opposite chiralit^EH. The 
simplest realization of a topological semimetal, one ex¬ 
hibiting only a single pair of nodes, necessarily breaks 
time-reversal symmetry. Indeed, T symmetry connects 
two Weyl nodes with the same monopole charg^^, im¬ 
plying t he ex istence of at least another pair with opposite 
chiralit}®^^. From this symmetry perspective, the coex¬ 
istence of 2D Dirac TI surface states and pairs of Fermi 
arcs is in principled allowed if T symmetry is respected, 
as inferred from ab-initio calculation^^. However, in the 
minimal two-Weyl-node model, this symmetry is broken 
and such coexistence seems to be mutually exclusive; the 
2D Dirac TI surface state is protected by T while Fermi 



FIG. 1. (Color online) Left: Phase diagram for the model § 
at 5o = 0 as a function of b = (5x,0,0) and M. The four 
gapped phases are a weak (WTI) and a strong (STI) topolog¬ 
ical insulator and trivial (I) insulator. The Weyl semimetal 
phases have six (WSMe), four (WSM 4 ^ 4 /) or two Weyl nodes 
(WSM2,2')- Solid lines represent gap closing or Weyl node 
annihilations at the corresponding Brillouin zone momenta. 
Upper right panel: Geometry, with periodic boundary condi¬ 
tions, used in numerical simulations. Lower right panel: A 
schematic of a proposed interference experiment to observe 
the coexistence of surface states. 


arcs are only realized in its absence. 

In this work we show how to circumvent this apparent 
dichotomy and realize both states at the same surface, 
the interface of a WSM-TI heterostructure, and we cal¬ 
culate the optical conductivity of the coexisting interface 
states, which serves as their distinct experimental signa¬ 
ture. We demonstrate this by modelling such interfaces 
using a canonical cubic lattice model describing TIs sup¬ 
plemented with symmetry brea king fi elds that can drive 
the system into a WSM phas^^^^. Spatially depen¬ 
dent parameter fields realize a generic model of a domain 
wall between two different phases. We numerically ob¬ 
serve that coexistence occurs when the Fermi arcs and 
2D massless Dirac surface states occupy distinct parts of 
the Brillouin zone delimited by the Weyl nodes^. This 
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observation is captured by a heuristic topological index 

Ja = Cana, ( 1 ) 

defined for each surface time reversal invariant momenta 
Aa and written in terms of known properties of the two 
phases (we use indices ... to label surface momenta 
and bulk momenta). Namely, the time reversal 

polarizations Tr^ = ±1 determine the p resence or absence 
of 2 D Dirac surface states at while Ca = 0(1) 

when a Fermi arc exists (is absent) in the vicinity of A^. 
By computing the index J'a for all one can predict 
whether coexistence of both surface states is allowed in 
a given surface or not (details are given below). 

Coexistence of Fermi arcs and 2D massless Dirac 
states —We model the bulk pha ses w ith the two orbital 
spinful cubic lattice HamiltoniarP^l^ 


H = ( 2 a) 

r zT 

+h-C. + My]4roCx, 

X,j ^ 

( 2 b) 

4 rWcx, ( 2 c) 

The position vector x runs over the sites of the cu¬ 
bic lattice, j = is a unit vector in each carte¬ 

sian direction, and fi = 0, x, 2 ). The operator Cx = 
(cxAt,CxA 4 ,CxBt,CxB 4 ) where Cxas annihilates an elec¬ 
tron in orbital a = A, 5 at x with spin s =t, The 
gamma matrices are given by T^ = {cFx ^SQ^dz^Sy^cTz^ 

Sx,CFy<S)So) and T^^^ = {cFy^Sz,(Jx<^Sx,CFx<^Sy,(Jo^Sz), 
where Sj and (jj are Pauli matrices describing the spin 
and orbital degree of freedom respectively and 5o, ctq are 
the corresponding identity operators. 

The first term H^i generally models a gapped insula¬ 
tor, apart from special parameter values that correspond 
to phase transitions between different insulating phases. 
At the time reversal symmetric momenta A^ of the cubic 
lattice, r = (0,0,0),X = 7^[(7r, 0, 0)], M = T^[(7r, tt, 0)] 
and R = (tt, tt, tt) with V the permutation operator, 
the gap is given by 2mi with mr,R = M ^ 3t and 
'^x,M = M ^ t. Depending on the relative signs of 
the masses one obtains a strong (STI, t < \M\ < 3t), a 
weak (WTI, |M| < t), or a trivial insulator (I, \M\ > 3t), 
cf. the horizontal axis in the phase diagram of Fig. 

The term Hb is parameterized by the four-vector by = 
(6o,b), where the pseudo-scalar bo breaks X and the 
pseudo-three-vector b breaks T. For now, we focus on 
the case of main interest bo = 0 and take h = bxX without 
loss of generality. We consider the effect of a nonzero bo 
later. With increasing bx the gap closes at one (or more) 
of the bulk time reversal symmetric momenta at which 
point the bulk spectrum is characterized by a 3D Dirac 
cone. Upon further increase of bx this Dirac cone splits 
into two 3D Weyl nodes and a WSM phase is obtained. 
Depending on the number of gap closings one obtains a 


WSM phase with two, four, or six Weyl nodes. A fully 
representative corner of the phase diagram is provided in 

Fig. [2 

To model an interface between two distinct phases we 
endow the parameters of the Hamiltonian ([^ with a po¬ 
sition dependence. We limit ourselves to sharp interfaces 
parallel to the x — z plan^^ with an infinite lattice in 
the X and ^ directions and a finite width Ly in the y 
direction. To avoid interfaces with the vacuum, we take 
periodic boundary conditions in the y direction resulting 
in two domain walls, see Fig. Explicitly, 

^^^f(6i,Mi) 

\{b 2 jM 2 ) otherwise. 

Depending on the values of ( 6 i,Mi) and ( 62 ,Af 2 ), we 
obtain an interface between any two phases of the model; 
In the following we focus on interfaces between a WSM 
and either a strong or a weak TI. 

In Fig. we plot cuts through the energy spectrum for 
kz = ^ and = tt as a function of kx for domain walls 
STI-WSM 4 and WTI-WSM 4 /, demonstrating the coexis¬ 
tence at zero energy of Fermi arcs and massless 2D Dirac 
states separated by the Weyl nodes. For the WTI and 
STI case there are an even and odd number of massless 
Dirac states at the interface respectively. The zero modes 
are doubly degenerate with one state localized at each 
domain wall. Importantly, not all STI-WSM or WTI- 
WSM domain walls have coexisting Fermi arc and Dirac 
states. For instance, at a STI-WSM 4 / interface (data not 
shown), where both types of surface states would have to 
exist in the same region of momentum space, only Fermi 
arcs are numerically observed. This suggests that coexis¬ 
tence at the interface is allowed only as long as the states 
do not overlap in momentum space. 

To assess the robustness of these states, we have stud¬ 
ied the effect of a finite X breaking term bo in the WSM 
side of the domain wall. For a bulk WSM nonzero bo 
acts as a relative energy shift for the Weyl nodes and 
may lead to the chir al magnetic eff ect, the existence of 
which is still debatecPinilltllEzEaHsni 

At the interface 

bo endows the Fermi arcs with finite dispersion and tilts 
the 2 D Dirac states, see right column of Fig. Increas¬ 
ing bo further ultimatel y ope ns up a gap in the WSM and 
destroys the Fermi arcs^^^. 

The coexistence is, furthermore, only allowed when b 
is perpendicular to the domain wall direction y; a finite 
parallel component {by 7 ^ 0) acts as a Zeeman term for 
the 2 D Dirac surface states and opens up a gap (as we 
have verified numerically). This effect is minimized by 
aligning the domain wall along b, which physically is an 
intrinsic magnetization and likely to be aligned with an 
experimentally identifiable crystallographic direction. 

Coexistence from bulk topology —Our numerical re¬ 
sults are captured by the topological index Ja defined 
in Eq. Q, the construction and use of which we now 
explain. The index relies on the observation that the 
bulk properties of the WSM and TI impose conditions 
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kz = 0.0, ho = 0.0 kz = TT, 6o = 0.0 kz = 0.0,6o = 0.25 



FIG. 2 . (Color online) Upper row: Band structure of a finite 
slab of WSM 4 -STI heterostructure with periodic boundary 
conditions [{bx,M)sTi = ( 0 , 2 . 6 ) and (&x,M)wsm 4 = ( 2 , 2 . 6 )]. 
Lower row: Band structure for a finite slab corresponding 
to a WSM 4 /-WTI configuration [{bx, M)sti = (0,0.2) and 
(bx, = (2,0.2)]. The = 0 , 7 r cuts shown in the 

first and second columns demonstrate the coexistence of 2D 
Dirac states with Fermi arcs in both cases. The third column 
shows the effect of bo = 0.25. Band structures are obtained 
for systems with linear dimension Ly = 160. 


on where their corresponding surface states must occur. 
When they are all compatible, the coexistence is allowed 
and protected by the presence of the Weyl nodes. 

To locate the TI Dirac surface states in momentum 
space we follow Refs. | 29 | and | 30 ] and define for each bulk 
time reversal invariant momentum the product of the 
parity eigenvalues 5 i of the filled Kramers pairs. For 
a surface perpendicular to a lattice vector, each time- 
reversal invariant surface momenta is a projection of 
two bulk momenta A^ and A^, see Fig.j^a). It has associ¬ 
ated with it the time reversal polarization tt^ = 5 i 5 j = ±1 
that determi nes th e number and connectivity of the Dirac 
surface state^^^^. Namely, for any path connecting A^ 
and Ka> with a ^ a' there is an odd (even) number of 
crossings at the Fermi level if TTa^Ta' = “1 (!)• For iLxi 
in Q ( 5 i = -sign[mi]. 

To similarly locate the Fermi arcs we employ the con¬ 
struction shown in Fig. [^b). Namely, we define a 2 D 
Hamiltonian by restricting the 3 D Hamiltonian ([^ to a 
cylinder enclosing a 3 D time reversal invariant momenta 
Ki but none of the Weyl nodes, and parameterize it by 
(ky^X) with A G [ 0 , 27 r] describing circles in the {kx^kz) 
plane. If a Fermi arc exists between the surface Brillouin 
zone projections of two given Weyl nodes around A^ it 
will cross this momenta as long as there is particle-hole 
symmetr^ F'^U^ I. For open boundary conditions in ky^ the 
2 D Hamiltonian has a midgap state at the intersection 
of the cylinder and the Fermi arc, see Fig. ib). This 
enables us to define an index Ca = 0(l^hat counts if 
there is a Fermi arc crossing A^ (or not).^ 


The two quantities Ca and Tr^ are separately obtained 
from the bulk of the WSM and the TI respectively. Since 
we are interested in a domain wall we combine them in 
the index Ja = 0, ±1 introduced in § associated to each 
interface momentum A^. From the index Ja we deduce 
the occurrence of nontrivial surface phenomena as fol¬ 
lows. At every A^ where Ja = 0 Si Fermi arc occurs. At 
a Aa for which Ja 0 the Fu-Kane criterion described 
above directly applies and analysis of the tt^ determines 
the existence of Dirac surface states. Hence, the key 
physical content captured by Ja is that both types of 
states at any given A^ are mutually exclusive. We note 
that Ca, and by extension J7a, relies on the fact that 
the Fermi arc crosses the 2 D cylinder shown on the left 
panel of Fig.j^b), which is guaranteed even if X breaking 
terms are present, as long as there is no gap opening in 
the WSM and particle hole symmetry is respected. 

To exemplify the use of 0 we apply the presented 
construction to the two domain walls considered above, 
as shown schematically in Fig. [^c). First, for the STI- 
WSM4 case we find that Ca = 0 for a = {kx, ky) = ( 0 , tt) 
and (tt, 0 ), which are thus intersected by Fermi arcs. Sec¬ 
ond, Ca = 1 and TTa = —1,1 for Aa = (0, 0), (tt, tt) re¬ 
spectively, indicating an odd number of crossings at the 
Fermi level between those two surface momenta repre¬ 
sented by a shaded circle in Fig. [^c) left panel. These 
conclusions are in perfect agreement with the numerical 
results shown in the upper panels of Fig. For the sec¬ 
ond domain wall, the WTI-WSM4/, the same procedure 
predicts two (an even number because of the WTI) Dirac 
surface states centered around (0,7r) and (tt^O) and two 
Fermi arcs crossing ( 0 , 0 ) and (7r,7r), see Fig. [^c) centre 
panel. Again, this agrees with the numerical results (see 
Fig. [flower panels). 

An immediate consequence of our analysis is that not 
all WSM-TI interfaces host coexisting surface states, 
even when b is aligned perpendicular to the domain 
wall. For instance, a domain wall involving the inter¬ 
polation (61, Ml) G STI ^ (^2,^2) ^ WSM4/ imposes, 
through Ja, that the Dirac nodes and the WSM Fermi 
arcs must cross = 0 at the same Aa, Fig. |^c) right 
panel. Since, as described above, Dirac states can occur 
only at Aa where there are no Fermi arcs, only the two 
Fermi arcs corresponding to the WSM4/ phase exist and 
cross {kx,kz) = (0,0), (7r,7r). Consistent with our analy¬ 
sis based on Ja, and as shown in the rightmost panels in 
Fig-i the inclusion of bo does not alter the coexistence 
of the surface states as long as it does not drive a phase 
transition to an insulator. 

Optical conductivity of the surface states —The pres¬ 
ence of the surface states alters the response to an ex¬ 
ternal electromagnetic field, which can be probed by op¬ 
tical spectroscopy. The reflection coefficients determine 
the optical response and are related to the optical con¬ 
ductive composed of a bulk and a surface state contri- 
butioiP^. Experimentally, the bulk optical signature of 
TIs has been accessed via optical spectroscopy^ and the 
bulk WSM optical signature is theoretically well under- 
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FIG. 3. (Color online) Construction of the topological in¬ 
dex: a) Projection of the different bulk band parity prod¬ 
ucts Si dehning the time reversal polarization TTa at each sur¬ 
face momenta Aa. The sign of al l Tig determine the pres¬ 
ence (absence) of TI surface state^^^^. b) A cylindric 2D 
Hamiltonian around a bulk time reversal invariant momenta 
Ai parametrized by ky and A G [0,27r) (left panel). For 
open boundary conditions in ky, the presence (absence) of 
a Fermi arc intersecting the cylinder, shows up as one (no) 
midgap state, defining Ca = 0(1) (right panel). Lower pan¬ 
els: Schematic interface Brillouin zone including the value of 
Eq. Q for three different domain wall conhgurations (from 
left to right): STI-WSM 4 , WTI-WSM 4 / and STI-WSM 4 /. 
The shaded regions enclose 2D surface massless Dirac states 
separated from the Fermi arcs (FA) in momentum space. 


stood^^^. Here, we compute the optical conductivity of 
the surface states reported above; in particular, the opti¬ 
cal conductivity of the Fermi arc, to our knowledge, has 
not been calculated before. 

We are interested in the linear response, long wave¬ 
length limit, where the incoming radiation has frequency 
uj and the momentum transfer satisfies p ^ 0. The Fermi 
arc and the massless 2D Dirac fermion exist in separated 
parts of the Brillouin zone. The total interface conduc¬ 
tivity is therefore given by the sum of their individual and 
independent contributions crsurf(^) ~ ^ ^fa(‘^)' 

The optical conductivity of a single 2D Dirac fermion, 
= cr|f (cd) = 1^, is well known from the con¬ 
text of graphen^^ and is isotropic and independent of 
the frequency uj. The Fermi arc can be modeled as a sin¬ 
gle chiral Fermion 7 /^+ with definite chirality, chosen to 
be positive without loss of generality. Its Lagrangian is 
C = To calculate its contribution we 

use the Kubo formula 

Re (w) = - lim -Im (w, p) (4) 

P^o UJ 

that expresses the optical conductivity in terms of the 
polarization tensor n*-^ (ci;,p). We find (se^^ for details 
of the calculation) 

RecT|f(a;) = Rea!J{uj) = Reagf(a;) = 0, (5) 


where 2k,q is the separation between Weyl nodes in 
momentum-energy space. This prefactor reflects the fact 
that the Fermi arc only exists on a bounded part of the 
2D surface Brillouin zone delimited by the surface Weyl 
node projection. The optical conductivity of the Fermi 
arc is thus highly anisotropic and divergent as cd ^ 0 in 
the clean limit. 


The total interface optical conductivity could 

be measured in a setup such as the one schematically 
shown in the lower left panel of Fig. where a TI thin 
film is deposited on top of a WSM, inspired by existing 
optical probeJ^. The TI bulk response vanishes for fre¬ 
quencies less tha n the bulk gap while the bulk WSM is 
proportional to coPUMl. Thus, for sufficiently low frequen¬ 
cies the response is determined by the surface and the 
coexistence can be probed by measuring the anisotropic 
Drude-like peak given in Eq. and a constant isotropic 
contribution from the 2D Dirac states. 


Discussion and conclusions —In this work, we have nu¬ 
merically demonstrated the possibility for Fermi arcs and 
2 D Dirac fermions to coexist at the same surface, the in¬ 
terface of a Weyl semimetal and a topological insulator, 
in spite of explicit T symmetry breaking. This is only 
possible if they do not coexist in the same region of re¬ 
ciprocal space and the time reversal breaking b-vector 
of the WSM is perpendicular to the domain wall direc¬ 
tion. We have introduced a heuristic topological index 
Ja , based on bulk topology, that can predict if and where 
the surface states are realized. This index captures the 
universal features of the bulk phases independent of the 
crystalline symmetry. Thus, it applies also to systems 
with rhombohedral symmetry, (e.g. the Bi 2 Se 3 family) 
that share the structure of the generic Hamiltonian stud¬ 
ied her^^. Even though Fermi arcs and Dirac cones can 
coexist, the latter are not as robust as those at Tl-trivial 
insulator interface. A component of the b-vector par¬ 
allel to the domain wall direction acts a Zeeman-term 
and gaps out the Dirac cone. The optical conductivity 
of the Fermi arc is found to be highly anisotropic and 
therefore optical spectroscopy serves as a probe of coex¬ 
istence of 2D Dirac and Fermi arc surface states. In sum, 
our results uncover the interplay of distinct topological 
bulk phenomena, topological insulators and semimetals, 
by showing that surface states with different nature can 
coexist at the same surface. 
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Appendix A: Appendix: Optical conductivity of the surface states 

In this appendix we derive the linear response of the surface states to an external electromagnetic field in the long 
wavelength limit, where the incoming radiation has frequency uo and the momentum transfer satisfies p ^ 0. Since 
the Fermi arc and the massless 2D Dirac fermion live in distinct regions of the surface Brillouin zone we calculate 
separately their contributions to the conductivity, labelled cr^^£(c(;,p) and <j£^(ci;,p) respectively. Mathematically 


j*(w,p) = = (T*']:(w,p) + (Tf^(w,p) £^j(w,p) 


(Al) 


where the current density j*(co’,p) and the external electric field F^^(cc;,p) depend on the external frequency uj and 
momentum p. 


For completeness and to establish notation, we begin by reviewing the calculation of the contribution from the 
single massless 2D Dirac fermion. This can be adapted from known results in the context of graphen^^ by simply 
dropping the valley and spin degeneracy factors. In linear response, we express the current in terms of the vector 
potential Aj as j'^ = Aj^ written in terms of the spatial component of the polarization tensor 11^^. To first order, 
this propagator is given by the bubble diagram shown in Fig. |^that represents the integral 


—G{k)rG{k + ph^, 


(A2) 


with G{k) = = (/co,k), and = (ci;,p). We use Feynman’s slashed notation defined through ^ The 

integral, written as 


n-(p) = f 

^ ’ J (27r)3 k^{k + p)^ 


(A3) 


has the precise same form as the QED(2+1) bubble and it is linearly divergent by naive power counting. Nonetheless, 
the final result is finite and can be obtained by using for example dimensional regularization. By analytic continuation 
to dimension d the result can be written in terms of the auxiliary Feynman parameter x G [0,1] as^ 




2e^ 


(47r)d/2 


Tr 


[ 1 ] /' 

Jo 


dx 


x{l — x)r{2 — d/2) 
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(A4) 


where Tr [1] is the trace over the dimension of the Pauli matrices. This last result will prove useful for the computation 
of the Fermi arc contribution to the conductivity. For the 2D massless Dirac cone one can set the spacetime dimension 
to d = 3 and Tr [1] = 2 to evaluate the integral and obtaiiP^ 
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(A5) 


We are interested in the real part of the optical conductivity, which can for example be probed by optical spectroscopy, 
related to the polarization tensor II*-^ through 


Recr*-^ (u;, p) = — lim — Imll*-^ (cj, p). 

P^O OJ 


(A6) 


Using that = 0 


Re (cj, p) = — lim -^Imfl^^fcj, p) = lim -^Im 

P^O PiPj p^O PiPj 


lb 


8 h 


(A7) 


where in the last step we have introduced the Kronecker delta function 6'^^ and restored e‘^/h to obtain the optical 
conductivity for 2D massless Dirac fermions used in the main texfP^. 

Now we want to calculate the contribution from the Fermi arc to the optical conductivity. The Fermi arc is a single 
chiral Fermion and thus its Lagrangian can be written as 


£ = + idi)'tp+ = 


(AS) 
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FIG. 4. Polarization bubble representing the integral (A2). 
evaluated at the corresponding momenta. 


The solid curved lines represent a fermionic propagator Gk 
For the 2D Dirac fermion it can be written as Gk = i/^ with = {ko,ki,k 2 ) 
while for the chiral Fermion describing the Fermi arc is given Gk = iV+lf^jk"^ with k — (A:o,A:i). The wavy lines 


represent the external electromagnetic field carrying external momentum p — (cj, p). 


where in the last step we have written the chiral fermion in terms of the projection of a Dirac fermion with the help 
of the projector operator V± = ^(1 ± 75 ). We note that the 7 matrices here do not correspond to the 7 matrices of 
the 2D calculation but instead can be chosen to have the representation 7 ^ = < 7 ^, 7 ^ = icjx and 7 ^ = 7 ^ 7 ^ = (Jz^ 

The propagator for such a state is Gk = with k = (ko^ki). Note that the chiral fermion is different 

from the usual Dirac fermion in 1 + 1 dimensions in that i) the propagator of the latter can be written as = 

iV-\-^lk^ + iV-^jk^ and ii) the former is embedded in a 2+1 dim ensio nal space time, unlike its 1 + 1 dimensional 
Dirac counter part. For the calculation of the polarization tensor (A3) of the Fermi arc, the integral in one of the 


directions is constrained by the separation between Weyl points in momentum space 2/^o, a fact taken into account 
by setting the limits of the first integral in k^ to —kq < kx < Therefore, for the Fermi arc case Eq. (A3) reads 




L ^ J w— k^kTw —• 

Using the definition of V-^ and the properties of the trace, the integral can be rewritten as: 
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(A9) 


(AlO) 


We can use the fact that 757 ^ = to write the second integral in the same mathe matic al form as the first by 

pulling the tensor out. Then, we notice that both integrals are formally the same as (A3) but in one dimension 
less. Thus, we can use (A4) with d = 2 and Tr[l] = 2 to write the final solution 


*nr;(p) = ^ - e^Ppppn] . 


(All) 


To relate this result with the actual surface conductivity notice that here /i, z/ G {0,1} while for the 2D Dirac cone we 
had G {0,1,2}. This means that the Fermi arc has only one non-vanishing longitudinal compon ent o f the optical 
conductivity (i.e. cr^^), while the 2 D Dirac fermion has both and contributions. Using (A 6 ), the optical 
conductivity of the Fermi arc is therefore 


e^Ko 


Re erg'(w) =0, 

ReafPco) = 0 . 


(A 12 a) 

(A 12 b) 

(A 12 c) 


The total surface conductivity is given given by the sum of (A7) and (A 12 a) for and by (A7) for a 





















